MACROSCOPIC DIFFUSION FROM A HAMILTON-LIKE DYNAMICS 
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Abstract. We introduce and analyze a model for the transport of particles or energy in 
extended lattice systems. The dynamics of the model has all the properties of Hamiltonian 
dynamics in a confined phase space : it is deterministic, periodic, reversible and conserva- 
tive. Randomness enters the model as a way to model ignorance about initial conditions 
and interactions between the components of the system. The orbits of the particles are 
non-intersecting random loops. We prove, as a weak law of large number, the validity of a 
diffusion equation for the macroscopic observables of interest for time arbitrary large, but 
small compared to the minimal recurrence time in the dynamics. 

3 

Fick's law of diffusion or Fourier's law of heat conduction describe phenomena part of 
everyday life : think of the diffusion of sugar in a cup of coffee or the exponential de- 
cay in time of the temperature of the same hot coffee. The laws of microscopic physics 
possess features that makes them look contradictory at first sight with the phenomenolog- 
ical laws of macroscopic physics. In particular, microscopic laws are reversible and when 
the dynamics is confined in phase space, the Poincare recurrence theorem ensures quasi- 
periodicity of the orbits. This a well-known problem which, in various guises, has gener- 
ated a large number of debates among physicists and mathematicians alike. Boltzmann 
made a decisive contribution to the issue by insisting on the fact that in large systems, 
the usual laws of macroscopic physics correspond only to a typical (with respect to the 
initial conditions) behaviour and not a uniform one. One should note that in general the 
argument can not rest on typicality of the initial conditions alone. There exist quite a few 
system for which macroscopic laws of normal diffusion are not obeyed : Fourier law is 
not observed in non-interacting gases, lattice of harmonic oscillators or Toda lattice. Still, 
given its pervasiveness, one can not expect "normal" macroscopic behaviour to depend on 
delicate microscopic dynamical properties. On the contrary, in the absence of any infor- 
mation on microscopic initial conditions or interactions reflected in a large scale structure, 
Fick's law or Fourier law seems to be the typical behaviour in Nature. Therefore, Boltz- 
mann solution to the problem should rather refer to typical initial conditions and typical 
dynamics. We propose an instructive toy model for Fick's law or Fourier law observed in 
systems which have a lattice spatial structure at the microscopic level . This model is de- 
terministic, reversible, periodic, conservative and amenable to a full rigorous treatment. 
Randomness enters the model as a way to model ignorance about microscopic initial con- 
ditions and microscopic interactions between the components of the system. With large 
probability with respect to this randomness, as the number of components increase, one 
can show that a Boltzmann equation holds and is equivalent to a macroscopic diffusion. 
The key point is that the minimal recurrence time in the dynamics increases as the num- 
ber of components does. This prevents any "recollision" occuring for time smaller than the 
minimal recurrence time. The model is reminiscent of the Kac ring model [1] but possess 
a rich orbit structure that makes it more similar to Hamiltonian dynamics. The orbits of 
the particles are non-intersecting random loops. 
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1. The dynamics and its first properties 
Le us consider a finite box of Z : 

A N = {i € Z : \i\ < N}. 

To each site of i 6 Ajy, we attach a ring 7£j carrying R sites fe 6 {1, . . . ,R}. The model 
consists of particles moving on 

C N = II U i = : ^ € * G A ^vi- 

The second ingredient of the model is the presence of "scatterers" that are located in 
between pairs of sites (k,i) and (k,i + 1). We define variables £,(k,i) taking values in 
{0, 1} and £(h, i) = 1 if and only if there is a scatterer between sites (k, i) and (k, i + 1). 
We put particles on the sites of Cat and a site carries at most one particle. At fixed intervals 
of times, all particles move forward on the rings. A particle located at site (k, i), namely 
at site k on the ring TZi will jump to site fc + lon ring 7^+i (resp. 7^_i), if and only if the 
following conditions are simultaneously satisfied. 

(1) There is a scatterer between (k, i) and (fc, (resp. (k,i — l)), namely £(k,i) = 1 
(resp. i — 1) = 1). 

(2) There are no other scatterer around that pair . 




Figure 1. A configuration of particles (black disks) and scatterers 
(crosses) on five rings. 

In every other cases, the particle located at site (k, i) simply moves forward on its own 
ring to (k + 1, i). Formally, we define the map r : Cn — >■ Cn giving the one-step evolution 
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of a particle located at (k, i) e Cn '• 

r(k,i) = J(k,i)(k + l,i + l) + J(k,i-l)(k + l,i-l) 

+ (1- J(k,i))(l - J(k,i-l))(k + l,i) (1.1) 

where 

j(k, i) = e(fc, - e(fc, i - i))(i - t(k, i + 1)) (i.2) 

and we set 

C(k,N-l)=C(k,N)=0. 

Here and in the following addition and substraction on the first component of a point in 
Cn are understood to be modulo R. The first factor in q 1 .2D accounts for the presence or 
absence of a scatterer, the two others are there to guarantee that r is an injective map and 
has thus orbits similar to Hamiltonian dynamics. 

Lemma 1.1. t is a well-defined bijective map o/Cn into Cn- t is therefore invertible. 

Proof. We first note that for any configuration £ of scatterers J(k,i) e {0,1}. Next only 
one of the three terms in (II . ID is non-zero. Indeed, assume that J(k, i) = 1 then the third 
term is obviously zero but so is also the second one because one must have i) = 1 
and thus J(k,i - 1) = £(& - l,i)(l - - 2))(1 - = °- Other cases are treated 

similarly. Thus, for any £, r(k,i) € {(fc+l,i+l), (fe+l,i — 1), (k+l,i)} C Cn and r is well- 
defined. We prove now that r is injective. Let us assume that we have x and x' such that 
t(x) = t(x') = (k, i), then both x and x' must belong to {(k—1, i), (k— 1, i— 1), (fc— 1, i+1)}. 
Suppose that x = (fc — 1, i), then because of the definition (I1.1D . we must have J(k — 1, i) = 
J(k — l,i — 1) = 0. This implies that r(k — l,i — 1) / (/c,i) and r(fc — 1, i + 1) / (fc, i) 
and thus that x = x' = (k — 1, i). Now, if x = (A; — 1, i — 1) and r(x) = (A;, i), then 
J(k - l,i - 1) = 1 and J(Jfe - l,t) = 0. Thus f(jfe - 1, t - 1) = 1 and {(k - l,i) = and 
therefore r(/c — 1, i) = (fc, i — 1) and r(A; — 1, i + 1) ^ (/c, i). Thus, in that case also, if x' 
is such that r(x') = r(x) then x' = x. The last possible case (when x = (k — 1, i + 1)) is 
treated in a similar way. We have thus proven that r : Cn — > Cn is an injective map. It is 
obvious that r is onto, we have thus shown that r is a bijection. □ 

One could add a degree of freedom ("velocity") to the model by deciding that all parti- 
cles move either in the positive (as above) or negative direction on the rings. We would 
then have two maps t + and t_ with r + = r and obviously r_ = t _1 . The dynamics is thus 
reversible in a analogous sense to Hamilton dynamics. 

To each x 6 Cn, we associate an orbit 

0{x) = {z € C N : 3n G N r n (x) = z} 

and a period 

T(x) = inf{n > : r n (x) = x}. 

Since the map r is injective, the orbits of two different sites are either identical or do not 
intersect. Moreover, any orbit is self-avoiding. 

Lemma 1.2. Every point ofCN is periodic (i.e. has a finite period) and we have, 

R < T{x) < R(2N + 1), VxsCat. 

Proof. Since the cardinality oi Cn is finite, the fact that r is injective implies that every 
orbit is periodic. This is actually a special case of the Poincare recurrence theorem and 
we do not repeat the proof. An upper bound on the period is given by the number of sites 
contained in Cn- A lower bound is obtained by observing that at each step one moves 
one step ahead on the rings. Thus, before coming back to the same site one must have 
performed at least one full rotation and this gives the lower bound on the period. □ 
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Another description of the dynamics is given by the occupation variable of each site 
(k, i) € Cn : a(k, i) 6 {0, 1}. The evolution of the occupation variables is given by 

a(k, i; t) = a(T-\k, i); 0), t € N* (1.3) 

or recursively by 

a(k,i;t) = (1 - J(k - - J(k - l,i - l))a(k - l,i;t - 1) 

+ J(k-l,i- l)a(k - 1, z - 1; t - 1) + J(/c - 1, i)a(A; - 1, % + 1; t - 1). 

(1.4) 

We note that at any time t e N*, the configuration a(-,t) is obtained as a permutation of 

the initial occupation variables <r(-; 0). Thus the dynamics is conservative. 

We sum up some obvious but crucial facts in the following lemma, for further reference. 

Lemma 1.3. (1) r~*(fc, i) depends on the configuration of scatterers only through the 
set of variables i) ■ £(k — n,i), n = 1, . . . , t, i = —N, . . . , iV}. 
(2) a(k, i; t) is a function of i) : £(k — n,i), n = 1, . . . ,t, i = —N, . . . , iV} and 

To summarize : for every configuration of scatterers, the dynamics has all the charac- 
teristic of a Hamiltonian dynamics in a confined domain : it is deterministic, reversible, 
(quasi-) periodic and conservative. 

2. Randomness and diffusion equation 

We let randomness enter the dynamics by taking a random configuration of scatterers 
and random initial occupation variables. Then r : Cn — > Ctv becomes a random map and 
(11 .3D implies that the configuration <r(-, t) at any time is simply a random permutation of 
a collection random variables. Likewise, orbits become non-intersecting random loops on 
Cn- Scatterers stand for interactions between individual components, namely the rings. 
In the absence of any macroscopic information on those interactions except its average 
"density", the most natural choice is to take the scatterers i) to be a set of independent 
Bernoulli variables of parameter < \i < 1. Similar considerations lead to take the initial 
configuration of occupation variables as a set of independent Bernoulli random variables 
but with a ring-dependent parameter to account for an initial inhomogeneous distribution 
of density of particles. We define the empirical density of particles on the rings IZi which 
are our main "macroscopic" quantities : 

1 R 

p R (i,t) := -J2°(k,i;t). (2.1) 
u k=i 

Our goal is to show that for times t smaller than the size R of the rings, this empirical 
density follows the solution of a diffusion equation, with large probability as R grows to 
infinity. As we indicate in the last section, it is easy to exhibit configuration of scatterers 
that would lead to an anomalous behaviour. 

Let p(i, t) be a solution of the system of discrete-time evolution equation : 

( p(i, t) - p(i, t - 1) = /x(l - /u) 2 (p(t - 1, t - 1) + p(i + l,t - 1) - 2p(i, t - 1)) , |t| < N 
) p(—N , t) - p(-N, t-l) = p(l- p){p(-N + l,t - 1) - p(-N, t - 1)) 
[ p(N, t) - p(N, t-l) = p(l- p)(p(N p(N, t-1)) 

(2.2) 

This is a discrete laplace equation in discrete time. Given initial data, the solution is 
unique for any positive integer t. We state now our main result. 
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Theorem 2.1. Let the {a(k,i,0) : (k,i) G C^} be a set of independent Bernoulli random 
variables and {pi : < pi < 1, i € A^} such that, E(cr(fe, i, 0)) = pi, Vfc G {1, . . . Let 
ako p(-,t) he the solution of the above system with initial condition p(i,0) = pi, Ve > and 
V0 < t < R a with < a < 1, 



lim ] 

R— too 



N 



(J {|p*(M)-p(M)|>e} 



-A' 



0. 



Proo/ By sub-additivity, it's enough to show 



lim F[{\p R (i,t)-p(i,t)\>e}] =0, ie Atv. 



(2.3) 



First, we show that 

E[p R (i, t)] = p(i, t), i G A N , < t < R a . (2.4) 

We consider the cases i / — N, N, it is easy to see that the cases i = —N, N are similar. 
Using (II .4D and noting that : 

J(k - l,i)J(k - l,i - 1) = 

and 

E[J(fc - 1, »)] = E[J(fc - 1, i - 1)] = p(l - p) 2 , Vl<k<R, 
we compute from (II .4D and (I1.2D : 

E[p fl (i, t)]-E[p i? (i, t-1)] = p(l-p) 2 {E[p R (i - 1, i - 1)] + E[p*(i + l,t — 1) — 2E[ / r R (i, f - 1)]) . 

We have used t < R a < R so that a(k,i;t) and £(k,j) are independent for all k, i,j by 
Lemma (11 .3D . Since E[p R (i,0)] = pi and because the solution to (I2.2D is unique, the 
relation (I2.4D is proven. We now look at the variance of p R (i, t). 



Var[p R (i,t)] 



fc=i fc=i / 

- / R R 

= R2 E i E ^ *)] - (E E [ a ( fc ' i; *)]) 2 1 

y fc,/c'=i fc=i 

We can write for the second term : 

E[ff(fc,i;t)] = E E k( T "*( fc ^);0)|r"'(A;,i) = x]P[r - *(M) = x] 

= ^ E[ff(a?;0)]P[T - *(M) = (2-6) 

Proceeding in the same way for the first one we obtain, 

E[a(k, i; t)a(k', i; t)] = E[cr(x; 0)er(x / ; 0)]P(r - * (Jfe, i) = x, r _t (fe', i) = x'). 



(2.7) 



When k ^ k',we get : 
E[<r(fc,i;*Mfc',*;*)] 



^ E[ct(x;0)]E[ct(x';0)]P(t" 



i) = x, r z) = x' 



(2.8) 
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because if k ^ k', then T~ t (k,i) / T~ t (k',i) for any configuration of scatterers, since r 
(and r -t ) is a bijection. The factorization of the expectation is obtained because the initial 
occupation variables are independent. Going back to (I2.5D . we take the first expectation : 

R R 
E[ E v(k,i;t)a(k',i;t)} = E E E[a 2 (x; 0)]P(t~* (fc, i) = x) 

k,k'=l k=lxeC N 

+ E E E[<T(x;0)]E[<T(x';0)]P(r-*(A;,i) = x, r"*^', i) = x') 
k^k'eC N x^x'&Cn 

(2.9) 

Since E[cr 2 (x; 0)] < 1, the first term above may be bounded by R and we get : 

R 

E E[a{k,i;t)a(k',i;t)} < E E[cr(x; 0)}E[a(x'; 0)]P[r-*(fc, i) = x, r _t (A;', i) = x'}. 

k,k'=l k+k' x^x'&Cn 

(2.10) 

For the second term of (I2.5D . we have : 

R 

(^E[crOM;t)]) 2 > E E E[o-(x;0)]E[<7(x / ;0)]P[T -t (jfe,t) =x]P[r -t (ife / ,t) =a/] 
fe=i fc^fc' x,x'eC N 

(2.11) 

Combining (I2.10D and (12.11ft , we get for the variance (1231) 

Var[p*(M)] < 1+ ^E E EK^;0)]EKx';0)]- 

fc^fc' x,x'GCjv 

(P[r-*(fc, i) = x, T~\k\ i) = x'\ - P[r _ *(A;, i) = x]P[r-' (A;', i) = x']) j 

(2.12) 

Using rotational invariance of both the distribution of the scatterers and the distribution 
of initial occupation variables, we get 



Var[p R (M)] < -+ -|E E Ek(*;0)]E[a(x';0)]. 

k'^lx,x'eC N 

(P[r - *(1, i) = x, t - *^, i) = x] - P[t - *(1, i) = x]P[t - *(A/, i) = x']) |. 

(2.13) 

Now, by Lemma (11.31) . ift+l<fe'<J?— 1 + 1 then T _i (0, i) and r~ t (k', i) are independent 
random variables and for those k', we have : 



P[r - *(1, t) = x, T~\k', i) = x'} - P[r"*(l, i) = x]F[T~\k' , i) = x'] = 0. (2.14) 
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We are thus left with 

R ' R 



Var[p R (M)]< -jU-^ £ E F[r- 4 (l^)=x,r-'(fc',,)=x'] 



Kfe'<t+1 



+ I E E P[r- t (l^)=^]P[r- t (fc',0=x'] 

fl-t+l<fc'<il X > X '€ C N 
l<fc'<t+l 

l + 4(*-l) 



6 

< —j^;, for i? large enough. (2.15) 

In the first inequality, we have used independance and E[<r(x;0)] < 1, in the second one 
the fact that probabilities sum up to 1 and the fact that there are 2t — 2 terms in each 
sum over k'. Finally, we have used the hypothesis that t < R a . We conclude by using 
Chebyshev's inequality and get : 

F[{\p R (i,t)-p(i,t)\>e}] <^U- (2-16) 
This in turn, implies (I2.3D . 

□ 

3. Relation to Boltzmann equation 

By using a "molecular chaos" hypothesis, one can derive heuristically a Boltzmann equa- 
tion for the evolution of the proportion of occupied sites on ring p R (i,t). This equation 
coincides with the diffusion equation (12.2ft . We start by an exact relation for the evolution 
of the number of particles contained in IZi, \i\ < N. 

R P R (i,t + 1) - Rp R (i,t) = X (i+1) _*(t) - X^ {i+1) (t) 

+ x (i_i)_n(t) - (t) 

(3.1) 

counts the number of pairs of sites on rings IZi and IZi+i for which, at time t, 
a particle will jump from ring IZi to 7£i+i and no particle jump from Hi+i to IZi. Namely, 

Xi_> i+1 (t) = \{k € {1, . . . , R} : a{k, i, t) = 1, a(k, i + 1, t) = 0, J(k, i) = 1}\. (3.2) 

We define also 

Xi_H+i(t) = \{ke{l,...,R} : a(k, i, t) = 1, a{k, i + 1, t) = 0}|. (3.3) 

All other X's are defined in an analogous way. In this context, the molecular chaos 
hypothesis amounts to assume : 

(1) The proportion X^i+iit) / R is independent of the scatterer distribution between 
the rings and therefore : 

1 R 

Xi^ i+1 (t)/R ~ ~Y,J(k,i) X^ l+1 (t)/R 

k=l 

~ E[J(k,i)\ Xi^ l+1 (t)/R 

~ /i(l-^) 2 l^+i(t) (3.4) 

as R oo. 
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(2) The proportion X^i + i (t)/R is given by the product of the proportions of occupied 
sites on IZi and vacant sites on TZi+i : 

X i ^ i+1 (t)/R = p R (i,t)(l-p R (i + l,t)) 

The first assumption is similar to the one made in the Kac model. The second one is more 
similar to a "genuine" Boltzmann assumption on the factorization of distribution of pairs 
of particles. One is thus lead to write : 

X^ {i+l) (t)/R = p(l - p) 2 p R (i,t)(l - p R {i + l,t)),i< N. 

A simple computation leads then to the evolution equation for \i\ < N. 

p R (i, t) - p R (i, t-l)=p(l- p) 2 (p R (i - 1, t - 1) + p R (i + 1, t - 1) - 2p R (i, t-1)). 

This is identical to (I2.2D for \i\ < N, equations for the evolution of the boundary densities 
are derived in a similar way. 

4. Conclusions 

We observe that taking a diffusive scaling limit (when N — > oo) of the discrete-time 
discrete-space diffusion equation (I2.2D will yield the usual diffusion equation 

d t p{x,t) = p{l - p) 2 d 2 p(x,t), x e [0,1], t € R+, (4.1) 

with Neuman boundary conditions. A more challenging task would be to derive 04. ID 
for a fixed parameter R directly in the diffusive scaling limit (when N — > oo). In that 
case, the absence of "collisions" of sufficiently distant paths expressed in (I2.14D is no more 
automatic and one should rely on a decay of spatio-temporal correlations. In order to 
obtain such a decay, it might be necessary to define the model on a higher-dimensional 
lattice. More generally, the possibilities of modifying the parameters of the model are 
numerous : one could also take correlated distribution of scatterers and study the effect 
of the correlations on the conduction properties of the model. 

We note that it is easy to find interactions (namely a fixed distribution of scatterers) 
between rings such that the empirical densities do not follow the diffusion equation (I2.2D . 
In contrast to the seemingly "chaotic" orbits that are typical of systems obeying normal 
diffusion, those configurations of scatterers tend to have a long-range order and give 
rise to very ordered orbits. In the set of possible interactions they represent "integrable" 
dynamics. By analogy, this might teach us a useful point to keep in mind when studying 
systems described by real Hamiltonian dynamics. We tend to conceive and study models 
that are (random or not) perturbations of dynamics that are "simple" or lend themselves 
to computations . But those systems might be rather untypical in the space of all possible 
interactions or indeed in real physical systems that obey the ordinary laws of macroscopic 
physics. 
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